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Abstract 

(N 

We demonstrate a new method for the calculation of inelastic scattering cross-section, which 
in contrary to the Regge-based methods takes into account the energy momentum conservation 
law. By virtue of this method it was shown that the main contribution to integral expressing 
inelastic scattering cross-sections comes not from the multi-Regge domain. In particular this 
leads to the fact that accounting of longitudinal momenta contribution to virtualities is sufficient 
and results in the new mechanism of cross-section growth. The necessity of taking into account 
I the large number of interference contributions is shown and the approximate method for this 

7-H purpose is developed. By considering the interference contributions from a single fitting constant 

^ achieved a qualitative agreement of the total and inelastic cross sections with experimental data. 

• I— I 
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1 Introduction 

Despite the fact the multi-peripheral model [2] has been used for description of hadron scattering 
for a long time, formal difficulties, which appear in calculating of inelastic scattering cross-section, 
in our opinion, are not overcame until now. These difficulties are caused by the fact that inelastic 
scattering cross-section with production of a given number of secondary particles in the finite state 

* E-mail: stiis@te.net.ua 
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(Fig.l) is described by the multidimensional integral of scattering amplitude squared modulus over 
the phase volume of finite state: 

x\T{n,puP2,...,Pn,Pl,P2,P3,P4)f6^^^ ^^3 + ^4 + Pfc " ^1 " ^^2 ^ (LI) 

where Mi and M2 are the masses of colliding particles with four-momentums Pi and P2; 
T{n,pi,p2, . . . ,pn, Pi, P2, P3, Pa) is scattering amplitude corresponding to inelastic process shown 
in Fig.l; 6^^^ is a four-dimensional delta function describing the conservation laws of energy and 
three momentum components in this process. Here it is also assumed that particles with four- 
momentums P3 and P4 are the same sorts as Pi and P2, respectively, and n secondary particles 
with four-momentums pi,p2, ■ ■ ■ ,Pn are identical. 

Since scattering amplitude is, in general, not a product of functions of some variables, and 
also due to the complexity of integration domain, the multidimensional integral in Eq.1.1 is not 
a product of smaller-dimensional ones. In considered inelastic process this domain of phase space 
of finite state particles is determined by the energy-momentum conservation law. As a result, the 
integration limits for one variable depend on the values of others. In order to overcome these 
difficulties one usually deals with the multi-Regge kinematics [5, 11, 8, 17, 4, 15, 10, 16]. The 
disadvantages of this approach we discuss in details in Section 7. 

The ultimate goal of this paper is to develop a new approach which is based on well-known 
Laplace's method [6] for the case, when scattering amplitude is set of multi-peripheral diagrams 
Fig. 1(b) within the framework of the perturbation theory. The essence of this method consists in 
finding the constrained maximum point of scattering amplitude squared modulus in Eq.1.1 under 
four conditions imposed by (J^^^-function of Eq.1.1. Then, expressing the scattering amplitude 
squared modulus as |rp = exp(ln{\T\'^)) , it is possible to expand the exponent of the exponential 
function in Taylor series about a point of the the constrained maximum, coming to nothing more 
than quadratic items. After that we obtain Gaussian integral, whose calculation is reduced to 
computation of matrix determinant of second derivatives with respect to /n(|Tp). 

The implementation of the aforementioned approach enables us to calculate the inelastic scat- 
tering cross section without the use of peculiar approximations of multi-Regge kinematics. As the 
result, we can draw the conclusion that the longitudinal momenta, which are usually neglected, 
play the significant role in the behavior of hadron inelastic scattering cross-section with energy ^/s 
growth. 

At the same time, according to the Wick's theorem, the scattering amplitude is the sum of 
diagrams of all possible orders of external lines attaching to the diagram in Fig. 1(b) (interference 
terms). In order to take into account these interference contributions one needs to modify the 
aforementioned procedure in the way, which will be outlined further in the paper. 

As will be shown, the value of inelastic scattering amplitude increases at the maximum point due 
to decrease of virtualities, which correspond to internal lines of "combs". Question, is this increase 
responsible for the total scattering cross-section growth, which is observed in the experiment, is 
also the subject of presented paper. 

The paper is structured as follows. In Section 2 we set ourselves the problem of finding the 
constrained maximum point of multi-peripheral scattering amplitude under the condition of energy- 
momentum conservation. Furthermore we consider some simplifications, based on the scattering 
amplitude symmetry properties. Afterwards, the analytical solution of the constrained maximum 
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(a) (b) 

Figure 1: 1(a) - a general view of an inelastic scattering diagram, 1(b) - an elementary inelastic 
scattering diagram in the multi-peripheral model ("the comb") 



problem is given in Section 3. The motivation of examination of the interference contributions at 
the calculation of inelastic scattering cross-section is discussed in Section 4. The calculation of 
inelastic scattering cross section, taking into account the contribution of all the interference terms, 
with the application of Laplace method at the relatively small number of final state particles is 
presented in Section 5. The approximate method for calculation of inelastic scattering cross section 
at any number of final state particles is developed in Section 6. Summary and conclusions are given 
in Section 7. 

2 The constrained maximum problem for the scattering ampHtude 
in multi-peripheral model. 

First, before we pass to the constrained maximization problem, let's examine the following simplifi- 
cations. According to Feynman diagram technique, the expression for scattering amplitude, which 
corresponds to a diagram in Fig. 1(b), has a form: 

T{n,P3,Pi,pi,p2,...,Pn,Pi,P2) = (-i<7(2vr)^)'(-iA(27r)^) 

xA{n,P3,Pi,pi,p2,...,Pn,Pi,P2) (2.1) 

with 

A{n,P3,P4,pi,P2r ■ ■ ,Pn,Pl,P2) = m^^iP^^Psy'-ie m^-{Pi-Pi.-piY -ie m^-(Pi^P^-p^-p2f-ie 

(2.2) 

1 1 

m^-iPi-Pa-Pi-P'i p„-i)^-ie rn^ - {Pi- P3-P1-P2 p„-i-p„Y -ie 

where 5 is a coupling constant in the outermost vertices of the diagram; A is a coupling constant 
in all other vertices; m is the mass of virtual particle field and also secondary particles. As in 
the original version of multi-peripheral model [2], pions are taken both as virtual and secondary 
particles. It was assumed that the particle masses with four-momentums Pi, P2, P3, P4 are equal, 
i.e.. Ml = M2 = M3 = M4 = M, where M is the proton mass. Note that the concrete choice of 
numerical value of mass M has no importance for the results presented in the paper. 

As it was noted in [4], for the most ratios of particle masses in the initial and final state, the 
virtual particles four-momentums on the diagram of Fig. 1(b) are space-like, i.e., their scalar squares 
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are negative in Minkowski space. The negativity of scalar squares of virtual four-momentums at 
the given mass configuration Mi = M2 = = M4 = M is easy to prove (see [13], page 30). 

Since the virtual particle squared four-momentums (Pi — Ps)^, (Pi — P3 — Pi)^, (Pi — P3 — Pi — 
P2 — ■ ■ ■ —pn-i —PnY' are negative at the physical values of four-momentums of final state particles, 
denominators in Eql.l are equal to zero nowhere in the physical region. Therefore it is possible to 
reduce ie to zero before all calculations. 

Due to negativity of virtual particle squared four-momentums the magnitude Eq.2.2 is real and 
positive. Therefore the search of the constrained maximum point of scattering amplitude squared 
modulus reduces to the search of the constrained maximum point of function Eq.2.2. Hereinafter 
we'll refer the expression Eq.2.2 as well as Eq.2.1, which differs from it by constant factor, as 
scattering amplitude, for short. 

Let us examine Eq.2.2 in c.m.s. of colliding particles Pi and P^. In such a frame of reference 
the initial and finite states have some symmetry, which is possible to use for solving the constrained 
maximum problem. In particular, the consideration of symmetries makes it possible to reduce the 
search of the constrained maximum of scattering amplitude to the search of the maximum of its 
restriction on a certain subset of physical process domain shown in Fig. 1(b). This restriction is the 
function of substantially smaller number of independent variables than the initial amplitude. 

For the further discussion of these symmetries and related simplifications, it would be convenient 
at first to take into account the conservation laws, expressing the scattering amplitude as a function 
of independent variables only. After decomposition of the three-dimensional particle momentums 
in c.m.s. frame to components, which are parallel and orthogonal kki. to collision axis, and lets 
name them longitudinal and transversal momentums, respectively. 

Energy of each particles in the finite state can be expressed by their momentum using the mass 
shell conditions, having n + 2 particles in finite state Fig. 1(b), that give us 3(n + 2) momentum 
components of these particles. Since we are looking for a constrained extremum, it is necessary 
to take into account four relations, which express an energy-momentum conservation law. It will 
result in the fact that amplitude Eq.2.2 can be represented as a function of 3n + 2 independent 
variables. The first 3n variables we choose are longitudinal and transverse components of momen- 
tums pi,P2, ■ ■ ■ ,Pn of particles produced along the"comb" in Fig. 1(b). The other two variables are 
the transverse components of momentum Ps±. 

If z-axis coincides with momentum direction Pi in c.m.s. and x and y axes are the coordinate 
axes in the plane of transverse momentums, the conservation laws look like 



-P30 + -P40 

^311 + ^^41! 
Pilx = — 
Pi±y = - 



= {PlO +P20 + ••• +Pno) 

= -{Pl\\+P2\\+ ■■■+P\\) 
{Pl±x + P2±x + ■■■+ Pn±x + P3±x) 
iPl±y + P2±y + ■■■+ Pn±y + ^4±j/) 



(2.3) 



where 



s = {Pi + P2)' 




(2.4) 



^M2 + (P3||)^ + (P3^,)V(P3^,)^ 
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Let's enter the following denotations: 



Ep = y/s~ (pio + P20 + • • • Pno) 
P\\p = Vl\\ +V2\\ + ■ ■ ■+Pn\\ 
Ppx = Pl±x + P2±x + • • • + Pn± 
Ppy = Pl±y + P2±y + • • • + Pn±y 



(2.5) 



X 



Then, solving the system Eq.2.3 for the unknown Ps\\, P4\\, Pi±x, P^iy gives: 



_ i;|-P||,^-P,\-2(p,,.P3x) / p , p L 4(AfH(P3x)^)(i^^P||/) 



(2.6) 



^tfk^z£Lz3i^k3A (p / 4(Af^+(P3x)^)(i^^P||,^) 

I ^l^^'^^V (E|-P||,^-P^,-2(P,,.P3.))^ 



where Pp± = {Ppx, Ppy) 

Substituting expression for P311 (Eq.2.6) into Eq.2.4 and resulting P30 into Eq.2.2 give us the 
scattering amplitude as a function of independent variables only, for which the conversation laws of 
all components of energy-momentum four-vector are taken into account. Below, referring to Eq.2.2, 
we assume that these substitutions have already been done. Taking into consideration this fact, we 
designate the amplitude Eq.2.2 as 

^ (n,p3±,Pi±,P2±, ■■■,Pn±,Pl\\,P2\\,--,Pn\\) (2-7) 

numerating only independent variables in the argument list. 

As it was shown in [13], for the further search of the constrained maximum we can limit ourselves 
to reduction of the scattering amplitude on a subset of the values of its independent arguments, 
which corresponds to zero values of transverse momentums of all particles in the final state. This 
reduction is a function of the longitudinal components of momentum Pi\i,P2\\, ■ ■ ■ ,Pn\\} which we 
designate as ^||(ra,pi|j,p2|j) • • • iPn\\)- Then from Eq.2.2 we get: 



(2. 



where pko = ^rin? + {pk\\Y, P30 = ^ + (^3,1)2. 

At the same time, assuming that all transverse momentums equal to zero, we have from Eq.2.6: 




- {Pio - Psof + {Pi\\ - Psi ^^ 

E Pko] + (Pl\\ - PsW - Y.Pk 
't=l J \ k=l 




^n-HP\\P±EpJl-^^JMl_^ 



^4|| = U P\\p^EpJ'^ 



(2.9) 



(^p)'-(^llp) 



Moreover, we have in c.m.s., Pio = -v/s/2) = — M^, where s is determined by Eq.2.4. 

For the further analysis it is convenient to switch from longitudinal momentums of secondary 
particles to rapidities yk defined by following relation: 

Pk\\ = m ■ sh ivk) , A; = l, 2,...,n (2.10) 
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n/2 particles 



(n-1 )/2 particles 



Symmetry axis of the diagram 



n/Z particles 



(n-1 )/2 particles 



Figure 2: An elementary inelastic scattering diagram in the multi-peripheral model with even (a) 
and with odd (b) number of particles on the "the comb" and it symmetry axis. 



Then function Ay can be written as = A^^{n,yi,y2, ■ ■ ■ ,yn)- The initial state in c.m.s is 
symmetric with respect to changes in positive direction of collision axis. In addition, those type of 
the diagrams presented in Fig. 1(b) have an axis of symmetry shown in Fig. 2 for the case of even 
number (left) and for the case of odd number (right) of secondary particles. 

It follows that for extremum search [13] we can examine the further reduction of scattering 
amplitude which is defined by the following equations 



^0 (n,yi,y2, ...,yR^ = Ay (^n, yi, 7/2, , -2/2, -yi) 



and at odd n 



^0 (n, 2/1,2/2, ...,ynzJ:) = ^11 ("-,2/1,2/2, ...,yzizj.+i = 0, -yn-i^i, -yn-i , -yi^ 



(2.11) 



(2.12) 



Considering the formula Eq.2.8 on subset, where reduction ^0 is considered, we will have Pp = 
by virtue of Eq.2.5. And therefore instead of Eq.2.9 we have the following expressions: 



En 



4M2 



(Ep) 



2 ' 




4M2 
(Epf 



or 




4M2 



\2 ' 




4M2 



(2.13) 



(2.14) 



Let us take into account that if we decompose all scalar square terms in denominators of Eq.2.2 
they will include the following difference F^ll - P311 = - M2 - P311 and negative Eq.2.14, 

chosen as the in the end will give us greater value in the denominator than in the choice of 
Eq.2.13. Therefore it is naturally to suppose that the main contribution to cross-section Eq.1.1 
gives the range of constrained maximum point determined by the scattering amplitude, where 
P3II and P411 are given by Eq.2.13, but not by Eq.2.14. Hence, considering the expression Eq.2.11 
and Eq.2.12 for the reduction of the scattering amplitude at zero transverse momentum region, and 
performing further transformations, we assume that the magnitude P^ and related with it quantity 

P30 = ^^Jm^~+~{P^)^ are expressed in terms of the longitudinal momenta of secondary particles by 
the relation Eq.2.13. 
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Applying symmetry relation and conversation of energy, it follows that on subset, on which the 
considered amplitude reduction Aq is defined, the energy corresponding to the line connecting n/2 
and n/2 + 1 vertices of the diagram in Fig. 1(b) is equal to zero in case of even number of particles 
at any values of independent variables (on which Aq depends). Similarly, for an odd number of 
particles the energy transferred along the line, which joins (n — l)/2 and (n — l)/2 + 1 verticies, is 
equal to m/2. The corresponding proof is given in [13]. 

Taking into account these results, the reduction of for the diagram in Fig. 1(b) with even 
number of particles can be written in the form, which is convenient for the further numerical and 
analytical calculations: 



/ — 



/ H 
' 2 



.2 ■ - • ■ 



X n I - E i^fc 1 +\sm-y. Pk\\ 

j=2 y \k=j J \ k=i 

where Sm = \/s/4 — — i-3||, Ek = m ■ ch (yk) and defined by Eq.2.10. 
The similar expression in case of odd number of particles in comb looks like: 




2> 



-1 



(2.15) 



^0 (n,yi,y2, ■■.,yn^ 




X m 



(2.16) 



X n - f + E ^fc + - E Pk\\ 

3=2 y \ k=j J \ k=i 

As it follows from Eqs.2.15, 2.16, it is convenient for the further calculations to make all quan- 
tities dimensionless by the mass m. In dimensionless form, these relations were used for numerical 
and analytical solution of the extremum for the reduction of the scattering amplitude ^o- The 
numerical solution of the constrained extremum problem for these expressions is described in detail 
in [13]. 

Summary results of the current section: It has been shown [13] that multi-peripheral scattering 
amplitude indeed has a point of constrained maximum under conditions imposed by the energy- 
momentum conservation law. The aforementioned symmetry relations for the constrained maxim 
point takes place. The rapidities of final state particles at the constrained maximum point produce 
an arithmetic progression. The difference of two adjacent rapidities on the "comb" increase with 
the growth of energy y/s. The main conclusion is that the absolute values of scalar squares of the 
final state particles four-momentums at the constrained maximum point decrease with the growth 
of energy ^/s, which leads to the growth of scattering amplitude at the constrained maximum point 
with the growth of energy. 

As will be shown further, this results in a new mechanism of inelastic scattering cross-section 
growth which has not been taken into account before. 
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3 Analytical solution of the constrained extremum problem as the 
approximation of equal-denominators 

We first consider in more detail the case of even number of particles n in the diagram of Fig. 1(b). 
The scattering amplitude reduction AQ(n,yi,y2, ■ ■ ■ ,yn/2) defined by Eq.2.15 undimensioned by 
mass m in this case looks like 



Ao(n,yi,y2,...,y|)= l-(Ec/i(yfc)j + (Pi3n)M U + ( AsH " (l/fc 



X n I 1 - f E c/i {yk)] + (Pi3\\ - E sh {yk) 

j=2 \ \k=j I \ k=l 




2N 

2 



-1 



(3.1) 



Here Piy = — and instead of a/s and M we use their dimensionless by mass m values. 

Since we search for the constrained extremum, under condition of energy-momentum conservation, 
it is assumed that P3 is described by Eq.2.13, in which again all values are undimensioned by mass 
m. In particular, taking into account the symmetry relation, normalization and the introduction 
of rapidity (see Eq.2.10), we obtain for Ep instead Eq.2.5: 

n 
2 

Ep = ^s-2^chiyk) (3.2) 
fe=i 

For further calculations we use the following denotation 

21 
2 

E = Y,ch{yk), (3.3) 

k=l 

AP = Pill - P311 (3.4) 

Note that due to Eqs.2.13,3.2 and Eq.3.4 quantity AP depends on rapidity as the composite 
function of E, which is denoted as AP{E) and the first term of Eq.3.1 depends on rapidity only 
via E. 

Instead of looking for the maximum of function AQ{n,yi,y2, ■ ■ ■ ,yn/2) '^e can look for the 
maximum of its logarithm, which we define as L: 



L = _2 In ( 1 - {Ef + {AP {E)f] - In | 1 + ( AP (ii;) - ^ iVk) 
-2 E In I 1 - f E c/i {yk)] + f AP {E) - ^ sh {yu) 

j=2 y \k=j j \ k=\ 

In addition, we make following denotations: 
Zi = \-{Ef + {AP{E)f 




(3.5) 



2 / . , \ 2 



Zj = \-\Y.ch{yk)\ + ( AP(ii;)-'E^^(2/fc) 1 ,i = l,2,...,f 

\k=j ) \ ^ k=i ) (3.6) 

%+i = l+ [AP(ii;)- t^s/i(yfc)J 
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Since after taking into account Eq.2.6, all variables of function AQ{n, yi,y2, ■ ■ ■ , yn/2) ^-nd hence 
logarithm became independent, then the extreme point can be found under condition that partial 
derivatives with respect to all variables are equal to zero. The equations for the extreme point 
problem can be written down in a form: 

'—1 " 
^ ^P{E)-Y.sh{yu) AP{E)-tsh{yk) 

^ = ^sh (yi) + Ach (y,) ^ ^ + 2ch {y,) = (3.7) 

dyi dE ^ Zj Zn+i 



j=2 ' 

1 ^ 

f AP{E)-Y: shivk) AP{E)- E shivk) 

+4chim) E ^ + 2ch{m) ^ 

j=i+i ' ^+ 



2 



(3. 



II - (n) + (^t) E + {y,) ^ (3-9) 

2 J — 2 2 



where 1 = 2, 3, ... ^ — 1. Eqs.3.7-3.9 form the system of equations for the extreme point 
search. An approximate solution of this system is the purpose of this section. The simplification of 
this system of equations can be attained in approximation, which we call "the equal-denominators 
approximation". Detailed justification for this approximation is given in [13]. Thus, for the further 
analysis of that system of equations at the maximum point we use approximation, in which all the 
denominators are equal to each other. We'll define their approximate common value as Z, i.e., 

Zj^ Z,j = l,2,--- + l (3.10) 

In this case the system of equations for the maximum point takes form 

fi-^||l£ I'M..)) . AP(.,-|:.M..) I ^0 ,3.n, 



j=2 \k=j J j=l+l \ k=l 



(3.12) 



k=l 




ch lyn] 

ch (yk) I + ; I AP (E) - J] s/i (yk) 1 = (3.13) 



sh (yn) 



k=l 
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From approximation Eq.3.10, in particular, we obtain Z„/2 ~ -^71/2+1 • Taking into account 
Eq.3.6 will result in identity: 

n/2 

AP {E) -Y^sh ivk) = — ^ (3.14) 

Inserting Eq.3.14 to the system of Eqs.3.11-3.13 we will get 

n / 21 \ 

ch (y^ ) I 1 , . N I Z dL ch (y^) , , 
r^+T^shivk) +^7^ + =0 3.15 



I 



j=2\k=j J j=l+lY^'^[y^) k=j J 2sh(ynjsh(yi) 



n / n 



j=2 \fc=j 



2{sh {yv_ 



Eqs.3.15-3.17 form the system of equations for the search of constrained maximum point of 
inelastic scattering amplitude in the approximation of equal denominators Eq.3.10. To compute 
this system we consider Eq.3.16 with I = f — 1: 



n -1 / n 



^^^i^ f + 2^K^t) ] + f i + 2 E f E (^^) 1 = (3-18) 

Subtracting Eq.3.17 from Eq.3.18 we have 

} i + (yt) - 2ch (yA - )'{^^ = (3.19) 

It follows from this equation that 



+ 2sh yik 

2sh{yr,.\ V^a 



2ch (yn)+ 

Taking into account that hyperbolic tangent is a monotonous function along the whole real axis, 
we have from Eq.3.20: 

= 3y| (3.21) 

Note, that this result agrees with the numerical results shown in Fig. 5 of Ref.[13]. Now prove by 
induction that 

Tl 

y|_fc = (2fc + l)yn, A: = l, 2,..., --1 (3.22) 
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Eq.3.22 is already proved for k = 1, since it coincides with Eq.3.21. Suppose that this equation 
is true for fc = 1, 2, . . . , n/2 — I — 1 (i.e., at y„/2-i, yn/2-2i ■ ■ ■ ■> Vi+i) ^^^d prove that it is true for 
k = n/2 — I (i.e., at yi). Subtracting Eq.3.17 from Eq.3.16 we obtain: 



2 E 4c/i(yfe) 



chjyi) 



2sh[ynjsh(yi) 







+ ^2 sh {vk) 



(3.23) 



Note, that sums ^ ch{yk) and ^ sh{yk) include only those yk with respect to / + ! < j < n/2, 

k=j k=j 



which are covered by the assumption of induction. Then, from Eq.3.22 we get ^ sh{yk). It makes 

k=j 

possible to calculate the sums included in Eq.3.23, which after transformations has a form: 



th{yi) = th 



n 



i] + i]y 



so that 



y, = (2(f-/) + l) yn 
yR^i = {21 + l)yii 



(3.24) 



(3.25) 



i.e., it coincides with proved Eq.3.22. 

Thus, for the diagrams with even number of particles we have shown that in the approximation 
of equal denominators, the rapidities produce an arithmetic progression at the maximum point and 
the ratios of all rapidities to the minimum rapidity produce the sequence of odd integers. 

To determine the values of rapidities, which constrainedly maximize the scattering amplitude, 
we also have to calculate yn/2i ™ which terms all rapidities are expressed. This can be done using 
the same equal-denominators approximation Eq.3.10 approach. 

In view of Eq.3.22 computing the sums in Eq.3.17 we have: 



(3.26) 



Z dL [in -l)y 

2dE^ „/ , / 

2{sh{y!L 



Now we can compute derivative dL/dE taking into account Eq.3.5 for the magnitude L: 



dL 
dE 



AE 



Zi fr^, z, 2 

\ 



AP{E)- E sh{yk) 

k=l 



\ 



Zi 



+1 



5AP {E) 



(3.27) 



Using the equal-denominators approximation Eq.3.10, we obtain: 



'dE 



AE 



( -1 - 

t AP{E)-Y.sh{yk) AP{E)-ish{yk) 

^P \E) sr^ fe=i i fc=i 

Zi ^ 7.. 



i=2 



Zn 



+1 



V 



9AP {E) 

'dE 



(3.28) 
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Figure 3: Comparison of the results of the numerical solution of Eq.3.33 (solid line) with the results 
of numerical maximization (circles) of the magnitude Ay(n, ^/s) at n/2 = 10 for the different energy 
ranges [13], GeV: 5 ^ 16000 (3(a)); 5 ^ 2000 (3(b)); 5 ^ 200 (3(c)). Here, it is taken into account 
that Ay(n, ^/s) = 2ysi. The good agreement of results shows the applicability of the approximation 
of equal denominators Eq.3.10 resulting in the Eq.3.32. 



After transformations with respect to Eq.3.14 we get: 



ZdL (nyn J sh (^(n + 1) yn_j gp^^^ (^) ^ 

Substituting Eq.3.29 to Eq.3.26, we reduce the resulting equation to the form 

/ N 5^311 (E) / N 

sh ((n + 1) y^j + ch ((n + 1) y^j = (3.30) 

Derivative dP^{E) / dE can be computed from Eq.2.13 with allowance for Eq.3.2 and Eq.3.4. 
Then the expression for ^311, which is made dimensionaless by the mass, can be written down as: 




P-i\\ = \\\^-E] -M2 (3.31) 

where it is assumed that ^/s and particle masses M at the ends of the "comb" are undimensioned 
by the mass m (pion mass was set to m and proton mass - to M). 

Taking the derivative of Eq.3.31 and substituting it into Eq.3.30, we obtain an equation, which 
after simple transformations looks like 

^-E = M-ch({n+l)yn\ (3.32) 



Note that the rapidity corresponding to momentum P311 is equal to (n + 1)^^/2 1 as it follows from 
Eq.3.31 and Eq.3.32. This expression would be obtained from Eq.3.22, if one accepts k = n/2 in it, 
i.e., arithmetic progression Eq.3.22 lengthens by one term. One might consider that rapidities of 
particles on the "comb" edges at the maximum point "continues" an arithmetic progression formed 
by the internal particles of the "comb" . This fact once again indicates the close relation between the 
equal denominators approximation Eq.3.10 and the arithmetic progression production by rapidities 
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Energy, Vs[GeV] Energy, Vs[GeV] 



(a) (b) 

Figure 4: Comparison of the results of the numerical solution of Eq.3.33(solid line) with results of 
numerical maximization (circles) of the magnitude Ay(n, -y/s) at ^ = 20 (4(a)) and ^ = 25 (4(b)). 
Designations are the same as in Fig. 3. The good agreement of results shows the applicability of 
the approximation of equal denominators Eq.3.10 resulting in the Eq.3.32. 



at the maximum point. In other words, the arithmetic progression production by rapidities at the 
maximum point is the consequence of equal-denominators approximation. 

It is also possible to verify this approximation in the following way. Taking into account E = 

n 
2 

J2 ch{yk) = sh{nyjk) /2sh{yri) (as it follows from Eqs.3.22,3.4) we have instead of Eq.3.32: 

fc=i ^ ^ 



M-c/i((n + l)yzi) (3.33) 



^ sh [ny 



2 2sh (yn 



This equation does not admit exact analytical solution, and further we will consider an approximate 
solution of it. However, we can verify permissibility of the approximations made above, which 
resulted in Eq.3.33, by the numerical solving of this equation at the different energies ^/s and by it 
comparison with the result of numerical determination of the maximum point. The results of such 
comparison are shown in Fig. 3 and Fig. 4. 

As seen from Figs. 3, 4, the "exact" numerical solution of Eq.3.33 practically does not differ from 
the results of numerical computation. This is the evidence of the fact that equal-denominators 
approximation Eq.3.22 is admissible approximation. Now let us consider the approximate analytical 
solution of Eq.3.33. Note that function sh{ny:i) /2sh{yii) in Eq.3.33 changes slowly at small values 
of yn/2 and can be approximately replaced by n/2 at yn/2 ~^ 0- In this approximation we obtain 
the following solution: 

ya = arccosh — — (3.34) 

^2 n + 1 \ 2M J ^ ^ 

Comparison of approximate solution of Eq.3.34 with results of numerical computation are pre- 
sented in Fig. 5. The Fig. 5 shows that Eq.3.34 gives a somewhat overestimated value in comparison 
with numerical computation. It is natural, since using approximation sh{nyik)/2sh{yn) ~ ^ in 
Eq.3.33, we underestimate the function sh{ny^) /2sh{y^) and, in that way, raise too high hyper- 
bolic cosine on the right-hand side of Eq.3.33. 

Nevertheless, as evident from Fig. 5, the absolute uncertainty of approximation Eq.3.34 does not 
increase with the energy growth, and because y^ itself increases the relative uncertainty decreases. 
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Figure 5: Comparison of the results of the approximate solution of Eq.3.33 (solid line) with the 
results of numerical maximization (circles) of the magnitude at ^ = 10 (5(a)), (5(d)); ^ = 20 
(5(b)), (5(e)); ^ = 25 (5(c)), (5(f)). The range of low energies close to the threshold branch point 
(in dimensionless form) is shown on (5(d)), (5(e)) and (5(f)). 



We can explain it also reasoning from Eq.3.33. Since sh{ny^ )/2sh{yii ) becomes small in comparison 
with Mch{{n + l)y^) at sufficiently high energies ^/s (and accordingly), the accuracy of the 
approximation for function sh{ny it) / 2sh{y n) has no importance. By neglecting sh{nyii) /2sh{yji) 
in Eq.3.33 in comparison with Mch{{n + l)y^) and by neglecting n in Eq.3.34 in comparison 
with y/s, we obtain the same result. It means that approximation Eq.3.34 ensures the "correct" 
asymptotic of value y^ at high ^/s. 

Let us note some features of Eq.3.34. Firstly, the approximate solution of Eq.3.34 undimensioned 
by mass m has a threshold branch point at {^/s = n + 2M). This means that difference of 
an arithmetic progression of the rapidity has same feature, which maximizing the amplitude of 
inelastic process. Contribution of the considered inelastic processes to an imaginary part of the 
elastic scattering amplitude after calculation with the help of Laplace's method [6] will be in some 
way expressed in terms of difference of the arithmetic progression Ay. Therefore it is possible 
to expect that mentioned threshold feature via difference of the arithmetic progression will be 
incorporated into imaginary part of the elastic scattering amplitude, which is required by unitarity 
condition. Note that Eq.3.34 has logarithmic asymptotic at the energies substantially exceeding the 
threshold value and at Ay = 2y!i ~ n~^, that coincides with the results of numerical computation 
(see [13]). 

For the following type of diagrams Fig. 1(b) with odd number of particles the whole procedure 
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is similar to described above for diagrams with even number of particles. At first we take derivative 
from the logarithm of amplitude restriction Eq.2.16 with respect to all rapidities yi,y2, ■ ■ ■ , Vn-i . 

2 

After that it is possible to use the equal-denominators approximation 

Zi ~ ... ~ Z n-i ~ .^ 71-1 . 1 = Z (3.35) 

2 2 

From the condition of equality Zn^_, ~ Zn-i we will obtain relation similar to Eq.3.14: 

2 2 

Ch ( gt/n-l J 

^111 - ^311 - E (y^) = ' \ (3-36) 

i=i 2sh [2y^ ) 

As it was done for Eq.3.20 using same recipe for the derivatives of logarithm of the scattering 
amplitude restriction with respect to and yn^_, and taking into account Eq.3.36 we have 

2 2 

= '^ynizl (3.37) 

2 2 

Then it can be shown by induction that 

y^-^k-i) ^ky,^, fc = 1, 2, - 1 (3.38) 

Thus all the rapidities, leading to the constrained maximum of the considered scattering ampli- 
tude restriction, can be expressed in terms of yri^. Repeating the calculations, which were made 

2 

in order to obtain Eqs.3.26-3.33, we have for this rapidity in equal-denominators approximation: 



"(f- '^^^'jl} =M-ch((^ + l]y._^] (3.39) 
^ 2sh 



2 

In an approximation similar to ones, which results in Eq.3.34, we get 



2 ( -Js — 

yn-i = — — -arccosh ] (3.40) 



^ n+1 \ 2M 

From Eq.3.40 it is evident that in case of odd number n the rapidity common difference, which 
constrainedly maximize the scattering amplitude, also has a threshold branch point. Note, that 
difference of an arithmetic progression is equal to yn^ in case of odd n and is equal to 2yn in case 

2 2 

of even n, i.e., as agreed, the approximation for the common difference of an arithmetic progression 
Ay{n, y/s) is expressed by the same formula both for even and odd n. 

The obtained analytical results enables to see how does the mechanism of virtuality reduction 
"work" with the energy growth. It's easy to show [13] that the scattering amplitude at the point 
of constrained maximum is expressed like 

^(0),n ^ (1 + ^))-2 ^ ^^^^ n))-("-') exp (c(V^, n)) (3.41) 



where 



a{Vs,n) = S 



'^^'"^=((!^^rjii"i) ' ^ (3.42) 

c(Vs,n) = 2(l-{n-l) (^s/M)"^ (^(^s/M)^ - l)) 



2_ 

hi 
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The a{^/s, n) and b{^/s, n) determine the characteristic value of virtuality at the maximum point 
of scattering amplitude and c(-y/s, n) determines the variation of virtuality along the "comb". In 
other words, the following estimate takes place 



1 

2 



2 / _ 1 ^ 2 



< 



,(^/s/M)"+i - 1, 



< 



(3.43) 



where 



is the absolute value of virtuality corresponding to j-th internal line on the " comb" 
in the point of constrained maximum. 

As it follows from Eq.3.43 energy included in same from that it useful to rewrite in this form 

{yfs/M) = exp (^^^ In (\/s/M)^ (3.44) 

It is obvious that the growth of exponent with energy ^/s is much weaker than the corresponding 

decrease with the growth of number of particles n. Thus, one can see that at not very small n 

1 

the value of (\/s/M)"+i ~ 1 even at high energies (\/s >> M). As the result, the difference 
of energy and longitudinal momentum squares is at least not negligible with respect to transverse 
momentum for each virtuality on the " comb" . This result comes in contradiction with the statement 
that virtulalities can be reduced to transverse momentum squares, which is usually claimed in the 

2 

standard approach [5, 11, 8, 17, 4, 15, 10, 16]. Taking into account the growth of {y/s/M)^+'^ 
with energy growth, we see that virtuality at the maximum point really decreases and the 

maximum value of amplitude grows with the growth of energy ^/s. Note also that at not very small 

1 

n the (\/s/M)"+i is close to unity at rather wide energy range which results in the much steeper 
growth than the one which is attained in Regge-based theories [2, 5] and described by factor of 
ln"'~'^ {^/s/M). Moreover, the higher n, the wider is the energy range. Thus the asymptotic 
behavior for different n is reached at different ^/s which enables to doubt the validity of the 
asymptotic formulas of multi-Regge kinematics. 



4 On the need of taking into account the interference terms at 
the calculation of inelastic scattering cross-section 

Before we turn to further calculations we need to make an important remark. Let's recall the 
formula Eq.1.1. According to Wick theorem, the scattering amplitude is the sum of diagrams with 
all possible orders of external lines attaching to the "comb" . In the terms of diagram technique it 
looks as follows. Plotting the multi-peripheral diagram of the scattering amplitude Fig. 1(b) at first 
we have adequate number of vertices with three lines going out of it and n lines corresponding to 
the secondary particles as it is shown in Fig. 6 (a). 

"Pairing" some lines Fig. 6(a) in order to obtain the "comb", we will get a situation shown 
in Fig. 6(b). The weighting coefficient appearing from this procedure, is included to the coupling 
constant. Finally we have to "pair" the appropriate lines of particles in the final state with the 
remaining unpaired internal lines in the diagram of Fig. 6(b). 

If we marked by ii - the external line, paired with the first vertex; 12 - the external line, paired 
with the second vertex and etc.; then ik is an external line, which is paired with k-th vertex, so 
every diagram will be characterized by sequence ii, ^2, ■ ■ ■ ,in- And in this case the total amplitude 
is expressed by the sum of n! terms, each of them corresponds to one of n! possible index sequences 
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(a) (b) 
Figure 6: Drawing the diagrams of the "comb" type 



and therefore the inelastic scattering cross-section can be written as 

_ ((27r)^)'g-*A2" r. dP3 dPi yj dpk r(4) /^p , p : f> „, _ D _ p A ^ 

/ X ( \ ^^'^^ 

X E A{n,P:i,Pi,p,^,p^^,...,p,^,P2,Pi)\ A{n,P3,P4,pj^,pj^,...,pj^,P2,Pi) \ 

\Pih.i2,....i„) J \P(jlj2,...J„) / 

Designation Yl means that a sum of all terms corresponding to all possible permutations 

P(il,i2,...,i„) 

of indices ii,i2, ■ ■ ■ ,in is considered. Moreover, since function A is real and positive, the sign of 
complex conjugation terms in Eq.4.1 can be dropped out. Writing this expression in the form 



((2^)^)Va^" P rfP3 



Pijl,j2,--,jn) 



(4.2) 



xA{n,P3,P4,Pn,Pi2,--,Pi„,P2,Pl) A{n,P3,P4,,Pj^,Pj2,...,Pj„,P2,Pl] 



The integration variables in each summand of considered sum can be designated so that the indices 
ii,i2, . . . ,in produce the original placing 1,2, At the same time the indices ji, j2, ■ ■ ■ , jn 

will produce some replacing and summation must be made over all permutations of these indices. 
Regarding to this, we can write down instead of Eq.4.2 



" 4:ViPiP2f-iMiM2f ■' 2P3o(27r)^ 2P4o(27r)^ 2pok{27rf 
X5(4) I^Pg + P4 + p^p,, - Pi - $ (n,P3,P4,Pl,P2, ...,Pn,P2,Pl) 



(4.3) 



where 

^> (n,P3,P4,pi,p2, •••,Pn,-P2,-Pl) = ^(n,P3,P4,pi,P2,---,Pn,^2, A) X 

X E A{n,P3,Pi,Pj„Pj2,...,Pj^,P2,Pi) (4-4) 

-P(il.J2,---Jn) 

At Section. 5 we will compute the integer from Eq.4.3 apply the fact that the amplitudes A, 
entering in Eq.4.4, have the points of constrained maximum. 
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5 Computing the multi-peripheral diagram contributions to in- 
elastic scattering cross-section with the Laplace's method 

For the further analysis consider Eq.4.3 in c.m.s. framework, expanding three-dimensional particle 
momenta into longitudinal and transverse components with respect to the collision axis: 



^ ^ ((27r)-*) g-tA^" dP^^dP^^^ " '^PfcX'^Pfcll ^ 

x5 (^M^ + Pli + ^M2 + P2 + £ - ^/^) (5.1) 

X<J ( I] Pk\\ + -fsll + -P4II ) 5 ( I] Pfc±x + P3±x + ^'4±a; ) ( E Pk±y + -PsXj/ + -RiXy ) 
\k=l / \k=l / \k=l / 

The last three (5-functions, whose arguments are linear with respect to integration variables, we 
can vanish by the integration over P411 , P4±x, PiLy In order to take into account the rest 5- function, 
which expresses the energy conservation law let's replace ^311 by new integration variable 



\\M''+\Y.m +p-i\\j + \Y.PkL + p^Lj (5.2) 



Those, making the following replacement we must express P311 through Ep. The corresponding 
expression will coincide with Eq.2.6 with the positive sign in front of the square root. Moreover, 
introduce the rapidities instead of longitudinal momenta: 



Pk\\ = m± {pk±)sh{yk) 
m± {pk±) = \/m+plj_ 



(5.3) 



After these transformations we have 



^{n,yi,pi^,...,y„,p„^,Pm,P2\i,P3\\,P3±,P^,P4±) 



2^M^ + {P,^^)^ + {P,^f 



n X 

Ep=V^- Y, "i±kip±k)ch{yk) 

k = l 



(5.4) 



E mx(pfcx)sh(j/fe)+P3|| ],Pi±=-{ E Pk±+P3± 



where it is assumed that the magnitude P311 is expressed in terms of the other integration variables 
via Eq.2.6. 

Now turn to the dimensionless integration variables and made the following replacement pk± — t- 
P3_|_ — ). We designate the new dimensionless integration variables just as the old variables, 
for short. Moreover, replace expression for P311 by the same expression divided by m. The same 
concerns the constants in expressions for cross-section, i.e., the designations M and ^/s are used 
for proton mass and energy of colliding particles in c.m.s., which is made dimensionless with the 
pion mass m. 

Now introduce the following notations of integration variables in Eq.5.4 and designate the rapidi- 
ties yi, 2/2, . . •) yn as Xi, X2,. . ., Xn] x-components of transverse momenta of secondary particles 
Pi±x, P2±x,- ■ -jPnix as Xn+1, Xn+2, ■ ■ ■ , X2n', y-componeuts of transverse momenta of secondary 
Pi±y, P2±y,---, Pn±y as X2n+i,X2n+2,- ■ ■,X3n- Moreover, designate X^n+i as P3±x and X3„+2 as 

P3±y 
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As it was shown in the previous sections, that an integrand A{n, P3, P4,pi,P2, ■ ■ ■ ,P2i Pi,P2) in 
Eq.5.4, expressed as a function of independent integration variables, has a maximum point in the 
domain of integration. At the neighborhood of this maximum point it can be represented in the 
form 

. s / 3n+2 3n+2 . ^ / rn^\\ 

A{n,P3,P^,p^,P2,...,Pn,P2,Pi) = A("^^"{V^)exp(-^ D^b [X^ - xi°>j i^Xh - xl"^U (5.5) 

where A^'^^'^{y/s) is the value of the function Eq.2.2 at the point of constrained maximum; 
Dab = — Qx dXf, (^^(^)) " derivatives are taken at the constrained maximum point of scattering 
amplitude. In other words, the real and positive magnitude A determined by Eq.2.2 is represented 
as yl = exp(ln(j4)), and the power of the exponential function is expanded into the Taylor series 
in the neighborhood of the maximum point with an accuracy up to the second-order summands. 
The accuracy of the approximation Eq.5.5 can be numerically checked in the following way. The 
function A, defined by Eq.2.2 can be rewritten in the new notation, which were identified previously 

A = A(n,Xi,X2,...,X3„+2) (5.6) 
Gaussian approximation of expression Eq.5.5 for function Eq.5.6 denote as: 

( 3n+2 3n+2 \ 
"2 E E - ^a^) (^^ - (5-7) 

a=l 6=1 / 

Since the functions Eq.5.6 and Eq.5.7 depends on large number of variables, it is impossible to plot 
them, so we have to introduce a new functions 

Fg (a, b) = A (n, X^^\ . . . , X^, ) + a, , . . . , xfl,,xf^ + 6, . . . , xf^ (5.8) 



4^)- (a, 6) = A(^) (n, xf\ . . . , X^.^xf + a, X% X^^.X^ + 6, . . . , xf^,) (5.9) 

Three-dimensional curves of these functions Eqs.5.8, 5.9 can be easily plot the in the vicinity of 
the maximum point (i.e., at the neighborhood of zero of variables a and b). The typical examples 
of such curves are shown in Fig. 7 and Fig. 8, where it is easy to see that the approximation Eq.5.7 
works well in the wide energy range. The results similar to ones in Figs. 7, 8 were obtained at 
different values of ^/s, i, k and n. As one can see from Fig. 7 and Fig. 8, that true amplitude and 
its Gaussian approximation Eq.5.7 start differ visibly only in the parameter region, which makes a 
insignificant contribution to the integral. 

Now let us proceed with identification A{n, P3, i-4,Pi,P2, • • • ,P2, Pi-,P2) in Eq.4.4 and define all 
possible arrangements of indices 1, 2, . . . , n by P^^\P^'^\ . . . , p("'). The function of variables X^, 
k = l,2,...,3n-|-2, which corresponds to arrangement pC), define as Ap(i) (n, Xi, X3„+2)- It 
differs from the function Eq.5.6 just by renaming the arguments, and therefore it also has a point 
of constrained maximum under condition of the energy-momentum conservation. The value of this 
function at the constrained maximum point is equal to the value of function Eq.5.6, i.e., and it is 
equal to A'^^^'^''{y/s) according to the replacement made above. Thus, if xf\x^\ . . . , Xn^ are the 
values of variables Xi,X2, ■ ■ ■ , Xn at the maximum point, now same values Xf\xf\...,x'^'^ will 
be the values of the variables Xj^ , Xj^ , • • • , Xj^ at the maximum point. Analogously X^+i , -^^^+2 1 ■ ■ ■ ,X.^ 
are the values of variables Xn+j^ , Xn+j2 Xn+j„ at the maximum point, and X2n+i ' ^2n+2 > ■ • • ; ^3n 
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(a) 



(b) 



Figure 7: Functions F"^-'^'^(a, 6) (dashed line) and F-^ ^ {a,b) (solid line) at energy ^/s = 5 
GeV and n = 10. The general image (7(a)), and the zoomed image (7(b)) at the neighborhood 
of maximum point. Clear that in an area that makes the most significant contribution to the 
integral, the scattering amplitude does not differ from its Gaussian approximation Eq.5.7, which 
demonstrates the possibility of applying the Laplace method. 



Figure 8: Functions F"7^'^(a,6) (dashed line) and (o,^) (solid line) at energy ^/s = 900 

GeV and n = 10. The general image (8(a)), and the zoomed image (8(b)) at the neighborhood 
of maximum point. Clear that in an area that makes the most significant contribution to the 
integral, the scattering amplitude does not differ from its Gaussian approximation Eq.5.7, which 
demonstrates the possibility of applying the Laplace method. 




(a) 



(b) 
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- for X2n+ji , X2n+j2 j • • • ) ^2n+j„ ■ For short we label the index of variable, into which the variable 
a goes at given rearrangement, as P^'^\a) i.e., the variable Xa replaced by the variable Xp(i){a). 

If we denote the matrix of second derivatives of the logarithm of the function Ap(i) at the 
maximum point by Z)^'*' , we will get the following approximation for the function Ap(i) : 

Ap(o (n, Xi,X2, X3„+2) = iV^) 

X exp -2 2^ ^pn)t„^ p(l)(h^ ^Pm(a) - ^p(OfM " x^ 



6=1 

Taking into account that Eq.5.10 depends on variables Xp(;)(^) and Xp(i)^^^ just as a function A 
depends on variables Xa and Xi, and the second derivative is taken at the same values of arguments, 
we have 

Using Eq.5.11 rewrite Eq.5.10 in more convenient form. For this purpose introduce the matrices 
p('), / = 1, 2, . . . , n! and by multiplying it with the column X of initial variables in Eq.5.7, we get 
a column in which the variables are arranged in that way so that in place of variable Xa became a 
variable Xp(i)^^y At next iteration taking into account Eq.5.11 one can rewrite Eq.5.10 in a matrix 
form in the following way: 

Ap(0 (n, Xi,X2, ...,Xsn+2) = (^/^) 

xexp(^-i (^l^(p«)^DP«X-2(lW)''i^pWl+ ^^-^^^ 

where is a column whose elements are the numbers X^\ a = 1, 2, . . . , 3n + 2 in the initial 
arrangement. Eq.4.4 can be rewritten in the form: 



(5.13) 



$ (n, Xi,X2, X3„+2) = (.4(0)." ( V^))' exp (x^) ^^X^ j 
X ^f] exp (^-iX^^Wx + (x(0))^y«x) 



= z)+ (p(0'\^2)pW (5.14) 



where 



VW = D + Dp(^) (5.15) 

If now we turn to Eq.5.4, we can see that all the other coefficients (except ^) under the integral 
don't change the values under the permutation of arguments. We replace these expressions by 
their values at the maximum point and take them out from integral. From this, we introduce the 
following notation: 



dEp 



(5.16) 

1 — 1 ^ ^ 



where p!^'!^ is the value of expression Eq.2.6 corresponding to particle momenta, for which the 
scattering amplitude has maximum, i.e., at the X^^ and nondimensionalized by m. 



21 



The expression for cross-section in this case can be written in the form: 



/ / \T \ n\ 3n+2 / / \ T 

:jW'"(V5)expf-(x(0)j DX(^nZI U (iX, exp f -iX'^D^X + j 



(5.17) 



n 



As the value ^ sh{y^^) + in Eq.5.17 is the negative value of longitudinal component of 

k=i " 

momentum P^^i^ taken at the maximum point, it can be replaced by due to the symmetry 
properties that have been discussed above. 

Multidimensional integrals under the summation sign can be calculated by diagonalizing of 
quadratic form in the exponent of each of them. Such diagonalization can be numerical realized, 
for instance, by the Lagrange method. The large number of terms in Eq.5.17 is substantial com- 
putational difficulty, which we overcame only for the number of particles n < 8. 

To represent results of numerical computations, it is useful to decompose Eq.5.17 in the following 
way: 

-(XW) DX^'AJ2 lldXae^p(-lx^D(^^X+[x('^) (5.18) 



, , r. (^W.n(^))V(0),n(^)/(")(^) 



L = ^ — (5.20) 

2(27r)3 \mj 

Note, that here and in the following sections we will use the "prime" sign in ours notation 
to indicate that we use a dimensionless quantity that characterized the dependence of the cross- 
sections on energy, but not their absolute values. 

The Eq.5.19 differs from the inelastic scattering cross-section an{\^) only by the absence of 

factor (m)^ (2(2^ (m)^) ' "^liich is energy independent and its consideration allow us to 
trace the dependence of inelastic scattering cross-section on energy y/s (Fig. 9 and Fig. 10). 

From Figs. 9 it is obvious that derivatives of cross-sections with respect to energies along the real 
axis are equal to zero at points corresponding to the threshold energy of n particle production, i.e., 
though the threshold values of energy are branch points for the cross-sections, they have continuous 
first derivative along the real axis at the branch points. 

Function crg{y/s) monotone increases in the all considered energy range as it follows from Fig. 10. 
At the same time from Figs. 11 one can see that fp {\/s) has drop-down sections. Moreover, 
even on those sections, where f^\^/s),n = 2-^5 increase, corresponding o"^(\/s) decrease. It 
makes possible to conclude, that amplitude growth at maximum point (which is the consequence 
of virtuality reduction) is generally responsible for the growth of inelastic scattering cross-section. 
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Figure 9: The calculated values of a'j^{^/s), n = 0,1,...,8 in the range of threshold energies for 
1, 2, . . . , 8 particle productions. Via o"g(-y/s) was denoted one of the contributions from the diagram 
(shown on the right) to inelastic scattering cross-section. 
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Figure 10: The calculated values of a'^{y/s) for the energy range y/s = 3 95 GeV 
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Figure 11: Calculated values of f^\\/s) determined by Eq.5.18 for the energy range ^/s = 3 95 
GeV. Comparison of these results with Fig. 10 shows that the relation Eq.5.19 for the dependence 
of partial cross section on energy -y/s is the most significant dependence A^^^''^{^/s) due to the 
contributions of the longitudinal momentum: at low n f^p\^/s) is increasing, but the partial cross 
sections Fig. 10 decreasing due to the fact that growth of A^'^^''^{y/s) is not fast enough, and at large 
n, the presence of decreasing sections f^p\\fs) have not inhibit the growth of cross-sections at the 
expense of growth • 

As it evident from Fig. 9 and Fig. 10 for some values of energy Eq.5.19 has a positive energy 
derivative and for some values of energy Eq.5.19 has a negative energy derivative. This makes a 
question. If we form from them a quantities 

^"'(V^)=E^'^^'n(^^) (5-21) 

and 

8 

<y" {^s) = Y,L''a'n{^) (5.22) 

n=l 

where L is defined by Eq.5.20, is it possible to choose the "coupling constant" L so that the value of 
Eq.5.21 has a characteristic minimum for the total proton-proton scattering cross-section? Answer 
for this question is positive (see Fig. 12), i.e., the curves agree qualitatively at the close values of L. 
The energy range shown in Fig. 12 takes into account all the inelastic contributions. We find indeed 
very interesting result, that curves presented on Fig. 12 and on Fig. 13, where calculated values of 
Eqs.5.21- 5.22 are given at L = 5.57, qualitatively agree with experimental data [21, 23]. 

In this paper we have examined the simplest diagrams of (j)^ theory and we intend to compare 
the qualitative form of these cross-sections with experimental data, but do not claim quantitative 
agreement. It is possible to hope that the application of similar computation method to more 
complicated diagrams in more realistic models will lead to correct outcome. 

As known, within the framework of Reggeon theory the drop-down part of total cross-section is 
described by the reggeons exchanges with interception less than unity [20, 24]. The cuts concerned 
with multireggeon exchanges with participation of reggeons with intercept greater than unity are 
responsible for the cross-section growth after the reaching the minimum [5]. 
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Energy, vs[GeV] 

Figure 12: Calculated values of a'^{^/s) ed, L = 5.57 in the energy range -y/s = 3 25 GeV. In the 
presented energy range the cross section qualitatively agrees with that observed in the experiment 
Ref.[22, 23]. 




Energy, vs[GeV] Energy, vs[GeV] 



(a) (b) 

Figure 13: Calculated values of a'^{^/s) at L = 5.57 in the energy range: 13(a) -y/s = 1.89 25 
GeV (starting from the threshold of inelastic scattering); 13(b) ^/s = 3 25 GeV (where increase 
of total cross-section with energy growth is clearly visible). Presented curves qualitatively agrees 
with that observed in the experiment Ref.[22, 23]. 

As will be shown in Section. 6, the accounting of (T^(y^) at n > 8 will not change the behavior 
of function a'^{^/s) Eq.5.20. This means that within the framework of given model the summation 
of multi-peripheral diagrams, when we compute the imaginary part of elastic scattering amplitude, 
will not result in power dependence on energy because since this dependence is monotonic. This, in 
turn, means that the corresponding partial amplitude has no pole! And this obviously differs from 
the results of standard approach in calculations of multi-peripheral model and from the results of 
Reggeon theory (see f.ex. [2]). The detailed discussion of the differences of our approach from the 
standard one is given in Section. 7 
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6 The approximate method for taking into account the interfer- 
ence contributions at high multipUcity of final state particles 

Accounting of the diagrams, which differ by the order of external lines attachment to the "comb" 
(interference contributions), is important and can be seen from the following concerns. For the 
relatively small number of secondary particles (n < 8) we are able to calculate all the interfer- 
ence contributions in the direct way without any approximations (as it was done in the Section 
5). Furthermore, one can evaluate the fraction of the contribution from the diagram with initial 
arrangement of external lines to the sum of contributions of all permutations of external lines. 
Finally, how it can be noticeably out from Fig. 10 of [12] this fraction is small and therefore, we 
can't restrict ourselves to accounting of the diagram with the initial line arrangement only at the 
calculation of inelastic scattering cross-section, as it is usually done at the standard approach to 
multi-peripheral model [2, 5, 11, 1, 9]. 

Each interference contribution can be computed numerically. However due to the huge number 
of contributions and large number of secondary particles n the direct numerical calculation of the 
sum of interference terms is impossible. In this section we present a method which enables us to 
overcome this problem. 

Let us represent an expression for the partial cross-section as a sum of "cut" diagrams in Fig. 14. 
The notation used in this section are the same as before with the exception of 

Xsn+l = 5 {P3±x — Pilx) X3n+2 = 5 {P3±y — P4±y) (6.1) 

As it was shown in previous sections, if we carry out calculations in the c.m.s. of initial particles, 
the maximum is reached at the zeros values of the transverse momenta and secondary particle 
rapidities are close to numbers that generate an arithmetic progression. Denoting the difference of 
this progression through Ay{n, ^/s) and the rapidity of particle, to which the line attached to the 
k-th. vertex of diagram in Fig. 6 corresponds, through we have: 

(0) _ f n+l 



vr = l^^-fcj Ay(n,s),A; = l,2,--- ,n (6.2) 

The interference contribution corresponding to total "cut" diagram, which refers to the j-th 
summand in Fig. 14, is proportional to an integral of the product of functions Eq.5.5 and Eq.5.10 
over all variables. Denoting an interference summand corresponding to the permutation Pj through 
an{Pj) and calculating its Gaussian integral (at the same time, other multipliers besides the squared 
modulus of scattering amplitude in an integrand are approximately substituted for their values at 
the maximum point), we have 



< (p.) = exp ( -i ( (a1]°))^ DO) Axf ) ] ] (6.3) 



1 / / . 



-1 



where we use the following designations: AX)^' = -Pri(x(o)), Z)(j) = [L)-^ + Pfp-'^P, 

scattering amplitude at the constrained maximum point) It should be reminded that here and 
further we use the "prime" sign to define the dimensionless cross-section, as in this paper we are 



2^y/s/A- {Ep/2) ^{Ep/2f - M^^ (here A (^X(O)) = ^(o)'" is the value of 
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Figure 14: Representation of the partial cross-section as a sum of "cut" diagrams. The order of 
joining of lines with four-momenta pk from the left-hand side of the cut is as following: the line 
with pi is joined to the first vertex, the lines with p2 is joined to the second vertex, etc. The order 
of joining of lines from the right side of cut corresponds to one of the n! possible permutations of 
the set of numbers 1,2, ... ,n. Where Pj{k), k = 1,2, ... ,n denote the number into which a number 
k goes due to permutation Pj . An integration is performed over the four- momenta pk for all " cut 
lines" taking into account the energy-momentum conservation law and mass shell condition for each 
of Pk- 

interested in the qualitative behavior of cross-section values, and we do not claim quantitative 
agreement. 

The essence of our method is as follows. The maximum in the right part of cut diagram in 
Fig. 14 is attained at X = (^X^^^^ . In other words, a maximum of function, which is associated 
with the right-hand part of cut diagram, can be obtained from a maximum of function, which maps 
with the left-hand part of cut diagram, by the rearrangement of arguments. Then the value of each 
interference contribution is determined by the distance between points of maximum in the right- 
hand and left-hand part of cut diagram as well as by the relative position of these maximum points, 
since in different directions contributions to scattering amplitude fall off with distance from point 
of maximum, in general, with different rate, and also by the relative position of proper directions of 
the matrices D and Pj DPj. In other words, multiplying Gaussian functions corresponding to the 
right-hand and to the left-hand part of interference diagrams in Fig. 14 each time we will obtain as 
a result Gaussian function, which has the proper value at the maximum point (which we call the 
"height" of the maximum) and the proper multidimensional volume cutout by resulting Gaussian 
function from an integration domain (which we call the "width" of the maximum). 

We assume that summands in Fig. 14 are arranged in ascending order of the distance between 
the maximum points in the right-hand part and left-hand part of cut diagram (we denote this 
distance through r) so that "cut" diagram with the initial attachment of lines to the right-hand 
part of diagram corresponds to j = 1. In other words, the line of secondary particle with the four- 
momentum Pi is attached to the z-th top in the right-hand part of cut diagram in Fig. 14. As follows 
from Eq.6.3, the interference contribution exponentially decrease with the growth. However, in 
spite of this the interference contributions do not become negligible due to their huge number, which, 
as discussed below, are increases very rapidly with growth. The value of is proportional to the 
square of magnitude Ay(n, y^), which, as was noted above, is zero on the threshold of n particle 
production and slowly increases with distance from this threshold. Therefore, for each number n 
there is fairly wide range of energies close to the threshold, in which the sharpness of decrease of 
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the interference contributions with the increase is smah in the sense that it is less important 
factor than the increase in their number. At such energies, which we cah "low", the partial cross- 
section cj„ is determined by the sum of huge number of small interference contributions. When 
the magnitude Ay(n, ^/s) is increased with the further growth of energy y^, the decrease rate of 
interference contributions increases, while the growth rate of their number with the increase 
does not change with energy. At such energies, which we call "high", the main contribution to the 
partial cross-section is made by the relatively small number of interference terms corresponding to 
the small r^, which can be calculated by Eq.6.3. 

If we compose the n-dimensional vector (we denote it through y^^^) from the particle rapidities 
Eq.6.2, which constrainedly maximizes the function associated with the diagram with the initial 
arrangement of momenta in Fig. 6, vectors maximizing the functions with another momentum ar- 
rangement will differ from the initial vector only by the permutation of components, i.e., these 
vectors have the same length. Consider two such n-dimensional vectors, one of which corresponds 
to the initial arrangement, and another - to some permutation, then in the n-dimensional space it 
is possible to "stretch" a two-dimensional plane on them (as a set of their various linear combina- 
tions), where two-dimensional geometry takes place. Therefore, the distance r will be determined 
by cosine of an angle between the considered equal on length n-dimensional rapidity vectors in the 
two-dimensional plane, "stretched" on them. An angle corresponding to the Pj permutation we 
designate through 9j, < 6j < tt. 

Thus, each of the terms in the sum Fig. 14 can be uniquely match to its angle 9j. At the same 

time the variable z = cos{9) is more handy for consideration than an angle 9j. Using Eq.6.2, can 

be shown that the variable z can take discrete set of values: 

12 , , (n-l)n(n + l) , 

zi = i--, , ^^J ; / = 0,1,---,^ ^-^ ' (6.4) 

(n — Ij n (n + 1) o 

Note that although the relation Eq.6.2 for the rapidities of secondary particles is satisfied with 
high accuracy at the maximum point, it is still approximate. This means that those contributions, 
to which matched one and the same value of variable z in Eq.6.2, in fact, matched a slightly different 
from each other values of z. 

Consequently, to such contributions correspond a similar but unequal to each other distances 
between maximum points in a "cut" diagram. In addition, this distance, as was discussed above, is 
not a unique factor affecting to the value of interference contribution. Therefore, if each interference 
contribution is associated with the value of variable z by the approximation Eq.6.2, it appears, that 
the different values of interference contributions correspond to the one and the same value of zi 
(Fig.15). 

Thus, while each contribution is associated to some value of variable z in the approximation 
Eq.6.2, the value of contribution is not the unique function of z. However, the sum expressing the 
partial cross-section cr„ can be written in the following way 



(6.5) 



\ J 

where AA^; is the number of summands to which the value zj = zi corresponds in the approximation 
Eq.6.2. The average value of all interference contributions in Eq.6.5 is already the unique function 
of zi- Therefore, we introduce notation 
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Figure 15: The dependence of interference contributions on zi at ^/s = 1000 GeV: (a) n = 8, 
(b) n = 9. Here and in subsequent figures the interference contributions divided by the common 



3n+2 3n+2 



multipher exp { — S ^a'^ ^ab^b'^ ] are indicated on the y-axis. Obviously, that to the one 



a=l b=l 

value of zi corresponds a lot of different contributions, as well as that the average values of the 
logarithms of these contributions are placed approximately on a straight line (see below Eq.6.11 
and Fig. 16). 



where {a!^(zi)) is some function, whose form at "low" energies can be determined from the following 
considerations. 

For any multiplicity n when the values of parameter I in Eq.6.4 are small and when number of 
corresponding interference contributions is relatively small, we can directly calculate these elements 
and their sum. Denote the maximum value /, for which all interference contributions are calculated 
through Iq. In particular, in this paper we managed to calculate the interference contributions up 
to /o = 6. Partial cross-section can be written as 

{n-l)n(n + l) 

^; = 4'^+4'^= E ^n{pj)+ E ^NiK('i)) (6-7) 

Zj = Zi, /=/o + 1 

V J 



where a'n is the sum of contributions sufficient at "high" energies, and cr„ is the sum of contri- 
butions sufficient at "low" energies. Thus, the difficulties in the calculations of the huge number 
of interference contributions mainly relates to the range of "low" energies and can be reduced to 
the approximate calculation of {(t'^{zi)) and AA'^;. The remaining part of the section is dedicated 
to this task. 

As follows from Eq.6.3, the exponential factor exerts the most significant effect on the depen- 
dence of {a'^{zi)) on zi. Note that the expression (^AXj^^^ D^^'^AX^^^ entering into the exponent 

in Eq.6.3 depends only on those matrix D^^^ components, which are at the intersection of the first n 
rows and first n columns, since all column Axj^^ components starting with n + 1 are zero, because 
they are the particle momentum transverse components at the maximum point. If we denote a 
matrix composed of elements located at the intersection of the first n rows and first n columns 
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of the matrix Z)^-') through Dy'^ and a matrix, which is obtained from the matrix t) in analogy, 
through Dy, we have 

D\^^ = {Dy'+PjD-y'P^Y' (6.8) 



At rather low energies (see [14]) eigenvalues of the matrix Dy^ are approximately equal between 
themselves. In such an approximation we obtain 

/ / nn2 / \y^''^^Sp{by] \ 

(a; (zO> = [a (xW)) V (^/i) exp ^-^^ (1 - z,) {w {zi)) (6.9) 



with 

(^'» = ^ / ...... <"°* 

Assuming that the multiplier {w{zi)) is weakly dependent on zi , we obtain 

\y(^^^Sp(Dy) 

(o-n (^0) = (o-n (zio)) exp I ^^^^ {zi - Z^J | (6.11) 

where is the minimum value of zi for which all interference contributions can be numerically 
calculated. Therefore, we can also numerically calculate the magnitude {cr'nizi^)). From Fig. 16, 
where the values of In {{Cnizi))) obtained by numerical calculation over all interference contributions 
and obtained with Eq.6.11 are compared, it follows that such an approximation is acceptable at 
"low" energies. 

In order to roughly estimate the function {'w{zi)) we can replace it by the Taylor series expansion 
taking into account the contributions no higher than linear as it was shown in [14]. We will designate 
the expansion coefficientds through wo and wi . The values of wq and wi are found by the calculation 
of {w{zi)) at zi close to 1 and (—1). So, we have the following expression instead of Eq.6.11 

(zi)) = (o-n (zio)) (wo + wi{l- zi)) exp | — ^ (z/ - Z^J | (6.12) 

Now let us turn to approximate calculation of AA'^; and use a new variable 

(0) 

5.(0) ^ Vk (g_13) 

Ay(n,s)y(^^±iM!^ 

where y^^^ are determined by Eq.6.2, Y^^\ k = l,2,...,n are considered as the components of 
vector Y^^\ which, as it follows from Eq.6.13 is of unit length. 

Thus, the angle 6j between the vector y^^^ = (jj^\y^\ ■ ■ ■ , yli^^ and vector P^r^ (y^*^-*) obtained 
by the permutation of corresponding components is the same as the angle between the vector 
y(0) = (Y}°\Yi°\. . . ,yi°^) and (Y^^^y Moreover, as it follows from Eq.6.2 

= = -y'X ■■■J^' = -y':L,v k = l,2,..,n (6.14) 
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Figure 16: Comparing the values of ln{{an{zi))) obtained by a direct numerical calculation with 
consideration of all interference contributions (circles) and by approximation Eq.6.11 (solid line) at 
n = 8, Vs = 10 GeV 16(a); 16(b) n = 9, = 10 GeV; 16(c) n = 8, ^ = 100 GeV; 16(d) n = 9, 
^/s = 100 GeV. 
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(a) (b) 

Figure 17: 17(a) - a sphere S2 and figure (is shown by points). Basis in the four-dimensional 
space is chosen so that one of the basis vectors coincide with vector 64 = (^i^i^iDi and the 
three basis vectors of three-dimensional subspace, into which depicted sphere is embedded, are 
perpendicular to 64; 17(b) - shortest arcs joining the points of figure ^41: into the two "hexagonal" 
and one "tetragonal" regions. 



It follows from here that all vectors A ^ ( Y^^^ ) are orthogonal to vector e 



, ,^ o l/Vn, l/\/n, ...,l/\/n 

J \ J 1 V ^ ' 

y n components ^ 

Therefore, considering vectors Pj^ (^^^^^ as the elements of n-dimensional euclidean space, which 

we denote through En, then the ends of all vectors Pj^ ^y*^*^)^ are lie on the unit sphere embedded 
into the (n — l)-dimensional subspace of We denote this sphere through Sn-2 and shape formed 
by the set of points in which the ends of vectors P~^ (y^^^^ (j = 1, 2, . . . , n!) come, denote through 
Fn\- In particular, when n = 4 the sphere S2 and figure F^i graphically look like in Fig. 17. 

The consideration of geometrical properties of figure Fn\ in [14] makes it possible to conclude, 
that at an arbitrary n a sphere 5^-2 can be divided into the parts of equal area, each of which 
contains only one point of figure Fn\ (as it shown in Fig.7(b-c) of Ref.[14]). 

Let us introduce a multidimensional spherical coordinate system so that the end of vector y^"-* 
is the "north pole" of sphere Sn-2- Then the number of figure Fn\ points, to which the values of 
variable z = cos{9) in the interval [z, z + dz] correspond, is equal to 

dN {z,dz) = p{z)dz (6.15) 

where 

Pi^) = ^f\^i^-^')^ (6-16) 

and r is the Euler gamma function. 

The fact that Eqs.6.15,6.16 are the appropriate approximations has been verified in [14] (see 
Figs.9, 10, 11 in [14]) 

Finally, using Eq.6.7, we are able to calculate the partial cross-section (T^(-y/s) at any number of 
final state particles n, where (cr'nizi)) is evaluated using Eq.6.12 and AA'^; is taken from Eqs.6.15,6.16. 
The comparison of such an approximation with the partial cross-section, calculated directly over 
all interference contributions for rather low n is given in Fig. 18. 
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Figure 18: The partial cross-section dependence on energy ^/s calculated over all interference 
contributions (solid line) and by Eq.6.7 with the application of approximations Eqs.6. 12,6. 15,6. 16 
(red circles): 18 (a) - a^(V^); 18(b) - a'gi^); 18(c) - a[o{V^); 18(d) - a[^{^); 18(c) - a[^{^). This 
approximation is acceptable at least in the range of parameters in which they are can be verified. 
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From this, we proceeded to the consideration of expression for the total cross-section and for 
the inelastic cross-section 

^'^ iV~s) = E (^^) (6-17) 

n=0 

'^max 

(^^) = E (^) (6.18) 

n=l 

which within the framework of the examined (p^ model is an analogue of total inelastic scatter- 
ing cross-section. Here Umax is the maximum number of secondary particles allowed by energy- 
momentum conservation law and A is the dimensionless coupling constant, which we considered as 
an adjustable parameter. Since the calculation of cr^ up to n = rimax takes a long time, we limited 
the upper bound of summation by those values of n, over which there are negligible contributions 
known to be smaller than the experimental error of cross-section measurements. 

Fitting constant L chosen to achieve a qualitative agreement o"'^(-y/s) and a''^{^/s) with observed 
in proton-proton collisions [21, 23] dependences on ^/s. The result of such a fitting presented in 
Fig. 19 and it qualitatively agree with experimental data not only at the high energies that is usually 
accepted in the Regge based theories, but also near the threshold of two-particle production (the 
first minimum of the total cross section Fig. 19(b)). This is due to the fact that the proposed method 
of calculation does not require any approximations, based on the asymptotically large energies. This 
may indicate that the experimentally observed behavior of cross sections is determined not by high 
energy asymptotic of the scattering amplitude as it is assumed in the contemporary approaches 
[24, 17, 18, 16]. 

However, the quantitative agreement with the experimental results [21, 23] requires the appli- 
cation of more realistic model than the self-acting scalar 0^ field model. 

7 The comparison of the presented approach results with the 
Regge-based ones. Discussion and conclusions. 

In the calculation of inelastic scattering cross-section with formation of n secondary particles a 
problem of the account for restrictions to a phase space, which come out of the energy-momentum 
conservation law, occurs. The traditional approach to this problem is based on the assumption 
that the main contribution to integral comes from the phase space domain, where the rapidities are 
ordered in such a way that each particle on the comb has rapidity much higher than the preceding 
one [16, 17, 18, 19]. Afterwards, some other auxiliary assumptions are usually made (see f.ex. 
comments to Eq.l of [9], Eq.l6 of [17] and Eq.4-5 [19]). 

As the result of these approximations one gets the energy-momentum conservation law that 
imposes restriction only on rapidities of the first and last particle on the comb. At the same time 
the rapidity of each intermediate particle on the comb varies from the value of the rapidity of 
preceding particle to rapidity of succeeding one at the integration over the phase space. 

Moreover, the difference of energy and longitudinal momentum squares is assumed to be neg- 
ligibly small (see comments f.ex. to A2->2+n in [18] and comments to Eq.l6 of [17]) in comparison 
to transverse momentum squared for each virtual particle propagator in the integrand of the equa- 
tion for scattering cross-section. Therefore integrand no longer depends on rapidities of final-state 
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Figure 19: Theoretical dependences of the a'^ {y/s) 19(a) and a'^{^/s) 19(b) obtained for the energy 
range ^/s = 1 100 Gev at L = 5.51. First minimum for the total cross-section can be obtained 
only when we take into account contributions from the high multiplicities. Experimental data for 
the inelastic 19(c) and for the total 19(d) pp scattering cross-section Ref. [22, 23] presented for 
qualitative comparison with the prediction from our model. Note: data-points for the inelastic 
cross-section, obtained from the definition ainei = crtotal — (^elastic- 
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particles (see Eq.23 of [17]). Finally one gets an equation similar to Eq.l of [9]. Lets give it here: 



ya and ?/{, are the rapidities of initial particles. According to this equation, one can see that, for 
example, the point: 



lies within the integration domain as well as point yi = y2 = ■ ■ ■ = Un = Ub- At these points 
all rapidities of secondary particles are equal to rapidity of one of the initial particles, which in 
fact, roughly violate the energy-momentum conservation law. Moreover, considering these two 
points one can say that their total energy-momentum vectors are not equal. Thus, in addition to a 
violation of energy-momentum conservation law, one gets that the total energy-momentum vector 
of the system varies while integrating over the phase-space. Finally, since integrand doesn't depend 
on rapidities of secondary particles (due to the considered simplifications), one gets the vicinities of 
aforementioned points give the same contribution to the integral (cross-section) as the other points, 
where energy-momentum is preserved. 

The analogical trick which leads to the similar results is employed in [8] as well. Namely, a 
similar assumption, that the main contribution to integral comes from a peculiar domain in the 
phase space (Eq.20 in [8]), is applied on the way from Eq.ll to Eq.22 in [8], in a transformation of 
a Delta-function, which stands for the energy-momentum conservation law. 

Contrary, in our paper we rigorously account for the energy-momentum conservation law at the 
calculation of the constrained maximum point of scattering amplitude Eqs. 2.3-2.6. 

Note, that the problem of the account of energy-momentum conservation law at the description 
of hadron-hadron interactions has been already raised in [7]. However, in [7] the problem was 
considered in the context of energy-momentum sharing between Reggeons in multi-Regge processes. 
At the same time, the aforementioned arguments enable us to conclude that the Regge-based 
approaches itself are in the poor agreement with energy-momentum conservation law. 

Moreover, the claim that the integration domain, where the rapidities are ordered in such a way 
that the rapidity of each particle on the comb is much higher than the rapidity of preceding one, 
put the major contribution to the integral expressing the inelastic scattering cross-section is only a 
hypothesis, as well as the proposition that the contribution of longitudinal momentums to virtuality 
can be neglected. In this paper, without any auxiliary assumptions, we found which area puts the 
main contribution to the integral (namely, the vicinity of the constrained maximum point). In [13] 
the point of maximum of integrand is calculated numerically without any additional assumptions. 
These results lead to us a conclusion that the contribution of longitudinal momentums to virtuality 
is substantial. 

Indeed, according to [9, 4], if one neglects the longitudinal momentum in the integrand of 
inelastic scattering cross-section, then its dependence on the energy of initial particles in c.m.s. is 
determined only by the volume of longitudinal phase space (using the terminology of [4] ) . Therefore 
the only mechanism of inelastic cross-section growth with the energy at such an approximation is 
related to a growth of longitudinal phase space. 

However, the existence of point of constrained maximum of the integrand leads to a fact that 
not the whole volume of phase space is substantial for the integral, but only the vicinity of the 
point of maximum. But even more important is that the account of longitudinal contributions to 
virtualities results in the dependence of integrand on energy of initial particles in c.m.s. which 





yi = Va, y2 = yi = ya,---,yn = yn-i = . . . = yi = ya 
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manifests itself in the dependence of square of amplitude module {A"'{y/s) in Eqs.5.19, 6.3) in the 
point of constrained maximum on ^/s. This leads to a mechanism of cross-section growth which, 
by definition, cant be taken into account while neglecting longitudinal components and therefore 
can't arise in the calculations based on the Regge kinematics. 

One can see from Figs. 12, 13, 19 that the aforementioned mechanism of virtuality reduction 
can be responsible for the total cross-section growth, which is observed in the experiment. In our 
opinion this is the role of longitudinal momenta in high energy hadron scattering, mentioned it the 
title of the paper. 
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